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A persistent spin helix (PSH) is a robust helical spin-density pattern arising in disordered 2D 
electron gases with Rashba a and Dresselhaus d spin-orbit (SO) tuned couplings, i.e., a = ±/3. Here 
we investigate the emergence of a Persistent Skyrmion Lattice (PSL) resulting from the coherent 
superposition of PSHs along orthogonal directions - crossed PSHs - in wells with two occupied 
subbands v = 1,2. For realistic GaAs wells we show that the Rashba and Dresselhaus /?„ 
couplings can be simultaneously tuned to equal strengths but opposite signs, e.g., ai = Pi and 
02 = —^ 2 - In this regime and away from band anticrossings, our non-interacting electron gas sustains 
a topologically non-trivial skyrmion-lattice spin-density excitation, which inherits the robustness 
against spin-independent disorder and interactions from its underlying crossed PSHs. We find 
that the spin relaxation rate due to the interband SO coupling is comparable to that of the cubic 
Dresselhaus term as a mechanism of the PSL decay. Near anticrossings, the interband-induced 
spin mixing leads to unusual spin textures along the energy contours beyond those of the Rahsba- 
Dresselhaus bands. Our PSL opens up the unique possibility of observing topological phenomena, 
e.g., topological and skyrmion Hall effects, in ordinary GaAs wells with non-interacting electrons. 

PACS numbers: 71.70.Ej, 75.70.Tj, 72.25.Rb 


Topological spin textures in crystals arise in connec¬ 
tion with the electron-electron interaction. Skyrmions 
in the fractional quantum Hall regime UM. magnetic 
and multiferroic systems Q exemplify spin patterns char¬ 
acterized by topological invariants associated with the 
nontrivial winding of the spins. Non-topological helical 
spin patterns, e.g., spin-density waves in metals can 
also occur. When coupled to conduction electrons, the 
emergent electrodynamics of the non-trivial spin textures 
gives rise to fundamental phenomena, e.g., the topologi- 
cal and skyrmion Hall effects in chiral magnets ■ 

Here we show that non-interacting 2D electrons in two- 
subband quantum wells in with matched SO couplings 
of opposite signs ai = /3i > 0, a 2 = —/32 < 0, can 
sustain a Persistent Skyrmion Lattiee (PSL), Fig. [I] This 
should allow the observation of fundamental topological 
phenomena in ordinary (non-magnetic) GaAs wells. 

The formation of this skyrmion lattice can be easily 
understood for ballistic electrons (later on we include 
disorder). For a single-subband well with ai = /3i, 
the Rashba-Dresselhaus Hamiltoninan is effectively ID: 
Hrdi = ^aiUyPx/h = giLBO'yBy/2, i.e., an electron- 
momentum (pa;)-dependent Zeeman interaction with a 
unidirectional effective magnetic field By {y || [110], 
X II [110]). Here ay is the Pauli matrix and pB the Bohr 
magneton. The corresponding quantum evolution oper¬ 
ator is: UBDi{t) = e-^ 9 nBrryByt/ 2 h _ ^-iayQix/ 2 ^ where 
X = Pxt/m*, Qi = Am*a\lh?, and m* the electron mass. 
Hence a spin-up electron injected at a: = y = 0 precesses 
around this By held such that 
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Figure 1. (Golor online) (a) Energy dispersion for a GaAs 
double well with two subbands (no disorder) and (b) its poten¬ 
tial profile and wave functions, (c) Calculated SO couplings 
vs. Vg: intraband (interband) Rashba (y) and Dressel¬ 
haus Pv (F). The dot-dashed vertical line (orange) indicates 
the crossed PSH symmetry point ai = Pi and 02 = —P 2 - 
(d) Energy contours; the arrows pointing along the orthogo¬ 
nal axes X (pink) and y (green) define the subband SO fields 
within subband 1 and 2, respectively, (e) Persistent Skyrmion 
Lattice (PSL) pattern in the 2DEG. The size of the circles and 
arrows denote {sz) and {sx,y), respectively. Blue (dark gray) 
circles stand for spins up and red (light gray) for spins down. 
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This spinor leads to a spin-density wave in the first sub¬ 
band s^(r), with oc sin(( 5 ix), Sy = 0 , oc cos(( 5 ia;), 
and pitch Qi. This helical pattern also arises in the 
presence of spin-independent disorder and (time-reversal 
conserving) interactions, and is known as Persistent Spin 
Helix (PSH) [1,0 (see Ref. [10 for ‘ballistic PSHs’). Ko- 
ralek et al. first observed a PSH via transient spin grating 
spectroscopy [HI ; Walser et al. imaged PSHs using ti ne- 
resolved Kerr rotation microscopy 12|. A single PSH is, 
however, non-topological. 

By considering a second subband with a 2 = —/32 < 0 
(Fig. 1), we can generalize Eq. ([T|) so that a spin up 
electron injected into both subbands evolves to 


Urdi ® Urd2 




/ cos{Qix/2) \ 

0 


s\tl{Qix/2) 

1 

—y 

cos{Q2y/2) 



\ism{Q2y/2) J 


( 2 ) 


where UiiD 2 = Here the second subband gives 

rise to a PSH with pitch Q 2 = 4 m*a 2 //i^ and spin density 
oc 0 , Sy oc sin(Q 2 y), and s^ oc cos(( 522 /), orthogonal to 
that of the first subband. These crossed PSHs form the 
unconventional pattern oc s^(r) -|- s^(r) in Fig. 1(e): a 
Persistent Skyrmion Lattice (PSL), that shows regions 
of zero and max/min spin densities characterized by a 
topological invariant (skyrmion number). 

The PSL texture inherits the robustness of the crossed 
spin helices, which are protected by the underlying SU (2) 
symmetry (within each subband) in lowest order of the 


J 


cubic and interband SO interactions. More physically, 
this robustness follows from the partial cancellation of the 
linear-in-momentum Rashba and Dresselhaus SO terms 
for ai = Pi and 02 = ~/ 32 , which renders unidirectional 
SO fields within each subband [Fig. 1(d)], and underlies 
the emercence of spin-conserved quantities in the sys¬ 
tem [ 1 , 0 . As the electrons move, they undergo spin 
rotations about orthogonal effective magnetic fields thus 
forming the skyrmion pattern in Fig. 1(e). Note that our 
PSL is identical to the ‘spin crystal’ of Ref. [l0 . 

We have also derived analytical expressions for the PSL 
spin density [Fig. HKe)] in the presence of disorder both 
(i) quantum mechanically 0 and (ii) via diffusive equa¬ 
tions 0. [I3 - [l9j[ . We show that intersubband-induced spin 
relaxation limits the PSL lifetime similarly (in magni¬ 
tude) to the cubic Dresselhaus term; PSLs are then fea¬ 
sible. With no disorder, our energy dispersions feature 
two Dirac cones at k = 0 Fig. [Ija), an anticrossing with 
spin mixing [significant in InSb, Fig. [2][a)], and highly 
anisotropic four-branch Fermi contours Fig. [ijd). 

Model Hamiltonian .— We consider a quantum well 
with two subbands. The two lowest spin-degenerate 
eigensolutions are (r|k,z/,cr) = e'^^'’^<.pv{z)\(Jz), z/ = 1,2 
and 0-2 =t, 4 ': with energies where 

k is the in-plane electron wave vector and is the z^th 
confined well level. Here we generalize the usual single¬ 
subband Rashba-Dresselhaus Hamiltonian for this two- 
subband case, which reads to linear order in k [see Sup¬ 
plemental Material (SM), Sec. (LA), for details [ioj ]] 


n = 


oy.i{(Tykx r^xky^ Pii^cTyky cTxkx') 
r](^(Jykx ^xky'j Y'(^(Tyky ) 


where ax,y are the spin Pauli matrices, kx y the wave 
vector components along the x || [ 100 ] and y |[ [ 010 ] 
directions, and a,y, the Rashba and Dresselhaus in¬ 
trasubband couplings, respectively, for subbands zz = 1 , 2 . 
Note that E g. (151) accounts for SO-induced intersubband 
couplings lol. l29i v ia the parameters 77 (Rashba) and F 
(Dresselhaus) [30|. Note that H describes two usual 
Rashba-Dresselhaus systems [the 2x2 upper left (oi. 
Pi) and lower right ( 02 , P 2 ) blocks] coupled via the in¬ 
tersubband “off-diagonal blocks” (ry, F). The energy dis¬ 
persions of H display anticrossing near fee, e.g., Fig. 2(a) 
for InSb wells. Similar dispersions (not shown) hold for 
a GaAs. As we show in the SM [Sec. (H.C)], for typical 
electron densities the Fermi wave vectors are such that 
kp «kc for GaAs and kp ^ kc for InSb wells. Next we 
use Lowdin perturbation theory to decouple (in orders of 
k or the subband energy separation Ae) the two Rashba- 
Dresselhaus blocks in Eq. ©; this procedure is valid for 
k = kp « kc as we discuss in the SM [Sec (I)]. 


rji^dykx o'xky) -f T'i^o'yky (Txkx) 

^2^k^ 0^2{o'ykx OTxky) -\- P 2 {^rryky efxkx) 


For convenience, let us first rotate the axes around z 
such that X ^ X \\ [110], y —t y || [IlO] (z |[ [001]) 
and then perform a spin rotation R = with 

dz= To lowest order, we find the uncoupled blocks 

Utz = -k (-Oy -k Pl^)axky -k {a„ + Pl^)aykx, (4) 

in which m* is the effective mass, Pl^ = Pi — P^ — Pi'’, 
with i = x,y] Pi = "/{vlkllv) (“bare” linear Dresselhaus), 
Pi ~ 'YTTn,zl‘2, rzi/ the subband areal density and 7 = 
11.0 eVA the bulk Dresselhaus constant [sij and Pl'^ is 
a function of Ae = £2 — £ 1 , Rv, V and F [Eqs. (S12)- 
(S13), SM]. As we show in the SM [see discussion after 
Eq. (S26)], Pl^ « Pi] hence in what follows we take the 
couplings Pi'’ in Eq. O to be Pv = Pi — Pi. This shows 
that the intersubband couplings (ry, F) essentially do not 
alter the PSH condition within each subband, = ±/3,y. 

Equation (|3|) shows that our two-subband well can be 
described (to linear order in k) as two uncoupled “copies” 
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Figure 2. (Color online) (a) Energy dispersions -Ek,Ai,A 2 
(scaled by a factor of 10 for visibility) along || [ 100 ] (or 
kx = ky) for an InSb double well. The black solid lines corre¬ 
spond to the uncoupled (77 = F = 0 ) bands and cross at 
kc- For 77, F 7^ 0 these bands anticross (dashed lines). Away 
from kc, the coupled and uncoupled cases coincide. The label 
sets (1, 2, 3, 4) and (5, 6 , 7, 8 ) denote Fermi points along 
kx at E'p and Ep, respectively. Panels (b) and (c) show (ct^) 
and {(Jy), respectively, along k^. The solid and dashed lines 
correspond to the respective energy branches in (a). 


of the usual single-subband Rashba-Dresselhaus model 
Hj/, with renormalized parameters. Each copy has SU(2) 
symmetry at Ofi, = ±/3i, Q. Next we show that the unique 
matching ai = /3i > 0 and 02 = —fi 2 < 0 occurs in real¬ 
istic GaAs wells [Fig. 1(a) shows the energy dispersions 
this case, see Eq. (S27) in the SM]. 

SO couplings for GaAs wells. — We self-consistently 
solve both Schrodinger and Poisson’s equations to ob¬ 
tain the eigenfunctions <Pv{z) of a GaAs double well 
(similar results hold for a single wide well). Fig. [TJb); 
see SM, Sec. (II). From ipv{z) we calculate the SO cou¬ 
plings, Fig. He). While ai and have opposite signs 
and are very sensitive to the gate voltage Vg across the 
well, Pi and P2 {Pi ~ P2) are practically constant (29| . 
At 14 = —10 meV (dot-dashed vertical line), we find 
cti = Pi = —0.2 = P2 = 1-45 meV A (black dots), which 
enables crossed PSHs as we discuss below. 

Robust eigenspinors even with disorder. — More 
realistically, we now consider an arbitrary time-reversal 
symmetric potential V (r) (spin independent), e.g., due to 
non-magnetic disorder Q , in Eq. O: nf,,y = n. + V{r). 
For Q!i = /3i > 0 and a 2 = -P 2 < 0 [Fig. 1(c)], ad¬ 
mits eigenstates of the form {r) = fy 

) and 7 />[J'’'(r) = fy) for subband 1, and 

^/;|'“(r) = ta:) and 7/72“’(r) = (p{r)e~''^^y/‘^\ fx 

) for subband 2. Here Qi, = Am*ay/hP {v = 1, 2, Qi > 0, 


Q 2 < 0) and I ta:), I J-a:) (| ty), I iy)) are the eigenvec¬ 
tors of ax (ay). The “envelope function” (p(r) satisfies 
(—/i^V^/2?n*-I-17(r)) (/3(r) = — Si, + 2a1m*/h'^) ip{r) 

[33 ■ Because [Hdis^^y] = 0 = 0, the eigen- 

solutions of 'Hdis possess (i) robust spin states against 
non-magnetic scattering (the spin and orbital variables 
factorize in 'Hdis)^ ('i) definite SO-induced spin- 
rotation phases dependent only on the distance ‘traveled’ 
along X and y, respectively. From (i) and (ii) we construct 
next “skyrmion states”. 

Persistent Skyrmion Lattice: Quantum approach. — 
Let ipiir) = (p(r) ty) + iy)) /a/2 

and 7 /; 2 (r) = p{y) {e^^^y/'^\fx) + ix))/y/^, 

i.e., stationary spin up states at r = 0 for each 
subband v = 1,2 at the Fermi energy Ep. The cor¬ 
responding spin densities s'^(r) within each subband 
are: s^(r) = i 7 /i|(r)cr 7 / 7 i(r) oc sin(Qix)x -I- cos{Qix)z 
and s^(r) = ^ifl{r)cnj} 2 {i‘^) oc sm{Q 2 y)y + cos{Q 2 y)z, 
with (T = axX + ayjj -I- a^z. These are orthogonal 
PSHs. Considering now the stationary superposition 
'0(r) = [ 7/1 (r) 0 7 / 72 (i’)]/a/ 2 at Ep (this is feasible 0111 
as we discuss later on), we can calculate the spin density 
s(r) = 7/)f(r)(l 0 dcr)'7/7(r) ( in units of h) 

s{x,y) = i|:/7(r)|2|sin(Qia;)x + sin(Q27/)y + 

[cos((5ia7) + cos((52?/)]-2|- (5) 

Interestingly, s(r) = (sAr) + s^(r))/2 forms a ‘persis¬ 
tent skyrmion lattice’ [ll|. Fig. H®), arising from two 
orthogonal PSHs, along x (1st subband) and y (2nd sub¬ 
band). We assume (/^(r) ~ (“weak disorder”), i.e., 
|(^(r)|^ = 1 in Eq. ([5]), to obtain Fig. He). Our PSL inher¬ 
its the robustness from its constituent persistent spin he¬ 
lices. In analogy to Ref. 1^, we can define ii = s/|s| and 
show that the PSL is characterized by a skyrmion num¬ 
ber over its unit cell area S: ^ f^n ■ (dxA x dyiijdxdy. 
Next we corroborate the quantum results presented here 
via diffusive equations. 

Semiclassical approach. — Following Refs. 00, 
we solve a set of diffusive transport equations for the cou¬ 
pled dynamics of charge n'^{x, y, t) and spin sf y ^{x, y, t) 
densities in subbands 12 = 1,2, valid in the weak SO inter¬ 
action limit {aiPPvkpT/h « 1 , t being the momentum 
scattering time [see SM, Sec. (HI), for details]. 

At ai = Pi and <^2 = —P 2 (symmetry point), the 
Fourier components of the spin density sf{qx, Qy, t) obey 

sj(q,t) = AJ,+ (q)e<+('i)* + Al_{q)e<~^^^\ (6) 


where AJ j_(q), j = x, y, z, are amplitudes set by the ini¬ 
tial conditions, uj], lo], ^ = —Dq^ — Tf,±Cf,qx = oj\., 
^l,± = ^l,± = -Dq^ - Ty ± C^qy = and ^ = 
= —Dq^, with D = VpTf2 the diffusion constant, 
vp = hkp/m* the Fermi velocity, q^ = q'^ + q"^, = 

Aaikl^T/m*, Cy = —Aa 2 kpT/m*, Tf = Saikf.r/hf, and 
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Ty = ?>a\k\<T/I t?. Hence given an arbitrary initial spin 
density s(q, t = 0 ) [or equivalently s(r, t = 0 )], we can de¬ 
termine Sj(q, t) = sj(q, t) + s|(q, t) and the spin-density 
profile Sj{r,t) = J Sj(q,t)e~’‘^ ’^dq at any time t. Next 
we discuss how PSLs can be realized. 

Excitin g a PSL via transient spin grating — The 
setup in 11 | can be implemented with crossed lasers 
of wave vectors q“ = q'^x and q^ = q^y. This cre¬ 


ates orthogonal spin gratings with an initial spin density 
Sz{r, 0 ) oc cos {q%x) + cos [q^y) and S 2 :(r, 0 ) = Sy (r, 0 ) = 
0. In real space, the resulting 2 -component sj -b s^ reads 

s, (r,t) = i + + cos(g“a;) 

+ \ cos (q^y). 

( 7 ) 

Similarly we find for x and y, respectively, 

s, (r,t)= i sin(q“:r), ( 8 ) 

sy (r,t) = - i sin (qly). (9) 


The spin grating experiment of till uses one laser and 
finds two decay constants for Sz in a one-subband well 
[Eq. ([6])]. Here we have two subbands and two lasers; 
hence we Hnd four time constants in each subband, two 
of which are equal [Eq. ([T])]; see also SM, Secs. (HI.B), 
(HI.C). Equations ©-(HI) show that the spin density ex¬ 
citation decays to zero as t ^ oo for arbitrary and 
However, when the laser wave vectors are tuned to match 
the pitches of the crossed PSHs, i.e., q“ = Qix, (ai = 
/3i), and q^ = Q 2 y, (a 2 = -? 2 ), we have w+(Qi,0) = 
( 0 ,( 52 ) = 0 and hence (r, t —>• 00) = sin ((5ia;)/4, 
Sy(r, t— >-00) = sin((52y)/4, and (r, t —>• 00) = 

[cos ((5ia;) + cos ((52?/)]/4. This is the PSL within the 
diffusive approach (cf. the quantum result in Eq. ([5])). 

Self-forming PSL upon photo-excitation. — A single 
PSH evolves from a uniform photo-excited spin-polarized 
density (e.g., sAr 0) = 1) in one-subband wells as 
demonstrated in [l2l. By the same token, a PSL (crossed 
PSHs) will also emerge in this setting, provided that 
ai = (3i and a 2 = —? 2 - Essentially, all Eourier compo¬ 
nents of S 2 (r, 0) = 1 decay to zero ast ^ 00, except those 
with the two “magic” q’s: qi = Qix and q 2 = ( 52 ?/, thus 
leading to crossed helices or a PSL [see SM, Sec. (III.D)]. 

Detrimental effects to the PSL — Our PSL so far has 
an infinite lifetime. It is known that in single-subband 
wells the cubic Dresselhaus (neglected so far) limits the 
lifetime of persistent spin helices. In addition, our PSL 
arises in two-subband wells and inter-subband spin decay 
may be an issue. However, we show that the spin relax¬ 
ation rate due to the interband SO coupling (Elliott-Yafet 
type) [ 3 ^ is comparable to that of the cubic Dresselhaus 
(D’yakonov-Perel type) I,El l33[ in limiting the PSL life¬ 
time [see SM, Sec. (IV)]. Furthermore, deviations from 


the PSL condition ai = (and/or q ;2 = —? 2 ) such 
that -b 6 v, with \di,lcxj\ 1 {v = 1 , 2 ), induce 

spin scattering with (golden-rule) rates ~ Sf,, i.e., spin 
dephasing vanishes in linear order in 61 , Isl. [35 j. Hence 


(8l,|^ 

I- 

1 11, 

hi 


To mitigate the stringency of the “a = /3” condition at 
a unique value, we note that both a and (3 can be varied 
simultaneously for a single sample - while still keeping 
a = (3 - over a wide range of electron densities in single¬ 
subband GaAs wells as shown in Ref. [sil- This allows 
for helices with gate-tunable pitches and ultimately to 
skyrmion lattices with controllable lattice constants (pro¬ 


vided the hndings in 3l| hold for two-subband wells). 

Band anticrossing & spin texture in k space .— We 
now turn to the effects of the interband couplings ry and 
r on the energy spectrum of % [Eq. ([3]), no disorder]. The 
solid lines in Fig. ©a) show the bands for the uncou¬ 
pled case r = 77 = 0 [see Eq. (S27) in the SM 0 ]. Both 77 
and P couple these bands with distinct ‘spin’ and orbital 
quantum numbers. Here we focus on InSb wells for which 
the SO coupling is stronger as compared with GaAs. The 
new bands Ek,Ai,A 2 (^ 1,^2 = ± 1 ) for non-zero P and 77 
display anticrossings around kc [36|, Fig. ©a) (dashed 
lines), and a strong spin mixing in (ctj). Figs, ©b) and 
He) [37|, near the anticrossing. This follows from an 
interplay of intersubband couplings: when either one of 
them is null, no spin mixing occurs as only same-spin 
branches couple in this case [see Eq. (S56) and discus¬ 
sion following it in the SM]. For completeness, in the SM 
we present the spin textures along the constant-energy 
contours Ek.Ai.Aa = Ep, E'p [Fig. ©a)]. 

Novel topological phenomena in 2DEGs ?— Similarly 
to chiral magnets, we conjecture that a PSL formed on 

can possibly 
Hall effects in 

ordinary GaAs wells. These phenomena could arise from 
two mechanisms: (i) the Lorentz force from the emer¬ 
gent magnetic held due to injected electrons following 
the real-space spin texture and (ii) the induced emergent 
electric held (Faraday induction) arising from the time- 
dependent topological hux of the drifting PSL, which can 
drive the PSL lattice sideways Q. Point (ii) is more likely 
to occur in our system [i^. Recent experiments 38. 411 


top of an electricall y d rifting Fermi se a 138 
lead to topological [13j and skyrmion [3£ 


have successfully demonstrated the electrically-induced 
coherent propagation of helices in GaAs wells . Fur¬ 
ther theoretical work similar to that in Ref. [l^ is needed 
to fully describe the quantum transport properties of our 
PSL on top of a drifting Fermi sea, which can unveil 
topological phenomena in ordinary GaAs wells. 
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